Dissipationless Phonon Hall Viscosity 
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We study the acoustic phonon response of crystals hosting a gapped time-reversal symmetry 
breaking electronic state. The phonon effective action can in general acquire a dissipationless "Hall" 
viscosity, which is determined by the adiabatic Berry curvature of the electron wave function. This 
Hall viscosity endows the system with a characteristic frequency, lj^; for acoustic phonons of fre- 
quency ij, it shifts the phonon spectrum by an amount of order {oj/lj^)^ and it mixes the longitudinal 
and transverse acoustic phonons with a relative amplitude ratio of w/oj^ and with a phase shift of 
±7r/2, to lowest order in uj/ujv. We study several examples, including the integer quantum Hall 
states, the quantum anomalous Hall state in Hgi-yMn^Te quantum wells, and a mean-field model 
for px + ipy superconductors. We discuss situations in which the acoustic phonon response is di- 
rectly related to the gravitational response, for which striking predictions have been made. When 
the electron-phonon system is viewed as a whole, this provides an example where measurements of 
Goldstone modes may serve as a probe of adiabatic curvature of the wave function of the gapped 
sector of a system. 



I. INTRODUCTION 

One of the most important discoveries in condensed 
matter physics has been that there are distinct states 
of matter that are distinguished not by their patterns 
of symmetry breaking, but by their topological order. ^ 
Such topological states of matter (TSM) cannot be de- 
scribed by local order parameters, but can sometimes be 
characterized by quantized topological responses to ex- 
ternal fields. For example, the quantum Hall states^''^ — 
the first topological states discovered in nature — can be 
characterized by their quantized Hall conductance. The 
three-dimensional time-reversal invariant topological in- 
sulators discovered more recently*"^ can be characterized 
by a topological magneto-electric effect^'*. More gener- 
ically, topological insulators in arbitrary dimensions can 
be characterized by topological responses to electromag- 
netic fields^. However, many TSM cannot be character- 
ized by electromagnetic response. For example in topo- 
logical superconductors the charge conservation symme- 
try is effectively broken, and the electromagnetic field is 
screened. Thus more generic response properties need to 
be investigated in order to distinguish different topolog- 
ical states. ^""'^^ 

In this paper, we propose a response property, the 
"phonon Hall viscosity," for gapped states that break 
time-reversal symmetry. For a quantum liquid, the 
viscosity tensor rjijki is defined by the linear response 
Tij = rjijkiVki, with T^- the stress tensor, and Vij — 
^{diVj + djVi) the gradient of the velocity field Vi. Usu- 
ally, a finite viscosity indicates dissipation in the system, 
similar to a finite resistivity. However the viscosity can 
have a dissipationless component, associated with the 
part of rjijki that is anti-symmetric under exchange of 
the first and second pair of indices. This Hall viscos- 
ity can only exist in a system that breaks time-reversal 
symmetry and is analogous to the dissipationless Hall 
resistivity. The Hall viscosity has appeared in the hydro- 
dynamic theory of the A-phase of He-3^'^ and was studied 
for quantum Hall liquids by J. E. Avron et. al}'^'^^ . It 



has since been studied for various (2 -I- l)-d topological 
states including IQH states^^, the (2-|-l)-d Dirac model^^, 
fractional quantum Hall (FQH) states^^"^° and Px + ipy- 
paired topological superconductors^^. The Hall viscos- 
ity provides a probe of gapped time-reversal symmetry 
breaking states in the charge neutral channel, which does 
not require charge conservation and thus may be a suit- 
able response for TSC and more generic TSM that can- 
not be characterized by topological electromagnetic re- 
sponse. In particular, it was recently proposed that the 
Hall viscosity of a rotationally invariant system is related 
to the angular momentum carried by each quasi-particle 
of the system^^'^^, which is in turn proportional to the 
"topological shift" of the topological field theory on the 
sphere^^. However, the Hall viscosity is only defined for 
a liquid in continuum, since the stress tensor is a momen- 
tum current which is ill-defined if continuous translation 
symmetry is broken by the lattice. The discussions of 
Hall viscosity in the literature have been treating the 
electron system as a liquid without explicitly consider- 
ing the lattice effects. This approximation is in general 
problematic since the Hall viscosity intrinsically depends 
on a length scale, and there are two natural length scales 
in a gapped system: a length scale associated with the 
energy gap, and a different length scale associated with 
the electron density. In general, the Hall viscosity will 
depend on both of these length scales and will therefore 
depend on non-universal short range physics. Even if 
the results in the continuum approximation are applica- 
ble in some cases, we are facing another problem of how 
to observe the Hall viscosity in general since it is diffi- 
cult to measure the velocity and stress of the electron 
liquid. To distinguish the Hall viscosity defined in this 
traditional way with the phonon Hall viscosity that we 
study in this work, we refer to the Hall viscosity of the 
continuum electron liquid as the gravitational Hall vis- 
cosity since the viscosity tensor of the electron liquid can 
be considered as a response to an external deformation 
of the spatial metric Qij^'^- 

To solve these problems of the gravitational Hall vis- 
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cosity, we alternatively define the phonon Hall viscosity, 
which is the adiabatic response of the electron state to 
the deformation of the crystal, i.e., to acoustic phonons. 
Instead of the stress tensor which couples to the defor- 
mation of the spatial metric and is only well-defined in 
the continuum, the deformation of the crystal and the 
electron-phonon coupling are always well-defined. The 
linear response of the electron liquid to the crystal defor- 
mation leads to a correction to phonon dynamics. The 
strain field Uij = ■^{diUj + djUi) plays the role of the 
spatial metric Qij, with Uj the displacement field of the 
nuclei. The phonon Hall viscosity is defined by the linear 

response ^^^^ = kjkiUki+VijkiUki with H the electron 
Hamiltonian depending on the lattice strain Uij, Xijki the 
elastic moduli, and ijijki the phonon Hall viscosity. The 
phonon Hall viscosity is well-defined for all gapped states 
with a regular lattice background, and is physically ob- 
servable through phonon dynamics. 

In the rest of this paper, we will first give a general def- 
inition of the phonon Hall viscosity in Sec. II and show 
how it relates to the adiabatic Berry curvature of the 
many-body electron wave function. In Sec. HI we study 
the phonon Hall viscosity in several example systems such 
as quantum Hall states, quantum anomalous Hall states 
and {p + ip) superconductors. We make comparisons be- 
tween the phonon Hall viscosity and the gravitational 
Hall viscosity in the continuum limit, and in certain ex- 
amples, we analyze features of the Hall viscosity that are 
universal or depend only on the low-energy theory of the 
electronic system. Finally in Sec. IV we will discuss the 
physical consequences of the phonon Hall viscosity, give 
order of magnitude numerical estimates, and discuss ex- 
perimental prospects for the observation of phonon Hall 
viscosity. 



II. DEFINITION OF THE PHONON HALL 
VISCOSITY 

A. Effective action of acoustic phonons 

In this section we will define the phonon Hall viscosity 
based on a generic discussion of acoustic phonon dynam- 
ics. The dynamics of acoustic phonons can be described 
using the long- wavelength effective action in terms of the 
displacement fields u{r), which describe the displace- 
ment of an atom from its original location. For an in- 
sulator, dissipation can be ignored at frequencies below 
the energy gap and the long-wavelength elastic response 
is described by an effective action for the displacement 
fields u, which can be obtained by integrating out the 
electrons: e-'^'-f/M = J Vc^Vce-^^'"''''''*l When the 
electrons form a gapped state, the phonon effective action 
is, to lowest order, 



Seff = \ j d^xdt{pdtUjdtUj - XijkidiUjdkUi), 



(1) 



where Xijki = Xkuj = Xjiu are the elastic moduli; its 
symmetry under interchange of the first or second pair 
of indices follows from invariance of the energy under 
rigid rotations. For a gapless state, such as a metal or 
a magnet, the phonon effective action will not be a local 
theory in terms of the strain fields. In the presence of 
time-reversal symmetry breaking, "Hall" viscosity terms 
are allowed: 



6Sh = ^ / df^ xdtrjijkidiUjdkUi, 



(2) 



where rjijki = —"UkUj- For an inversion symmetric sys- 
tem, this is the only additional term that can be added, 
to cubic order in momenta. Anharmonic corrections to 
the acoustic phonon dynamics go like {du)^, so they are 
implying that the Hall viscosity term may be dis- 
tinctly measurable since it is lower order in momenta. 
Impurities in the crystal may be difficult to treat, but 
as we discuss in the Section IV B, their contributions are 
not sensitive to the sign of the time-reversal symmetry 
breaking of the electronic system and may be separated 
from Hall viscosity contributions. 

Up to total derivatives, (2) only depends on riijki+r]kjii- 

^^H = ^ J d'^xdf\i]±^idiUjdkUi + ri-.^^idi{ujdkUi)], (3) 

where r]f^^i = \ {riijki ± Vkju)- The boundary terms may 
have interesting consequences for the surface waves of a 
medium with such a Hall viscosity, but we will ignore 
them in this paper. In some cases it will be conceptually 
more clear to write the above action in terms of the strain 
tensor Mi j = ^{diUj+djUi) and the rotation tensor m,j = 

^Sh = ^ j d'^xdt[r]fjliUijUki + Tq^ki^ijihu 



+ '^VijkiUijrhki], 



(4) 



where ijff^^i = rj^g = -r^f^f^-, r]f^ 



AA 
kl 



'Ijikl 



and rj^j'^i = rj°^^i = -rjf-f^. can all be deduced from 
Vijkl- For an isotropic three-dimensional system, rjijki 
must vanish. In two dimensions, r]^^ always vanishes; 
for an isotropic 2D system, or one with 7r/4 rotation 
symmetry, rj^^ and rj^"^ each reduce to a single num- 
ber: r?^ ^ v'J.y = v'.y'yy''''' and r?^ ^ v'^ - 
In this case, 



xxxy 



n 



yyxy 



SSh = 2y d'^xdt[r]"{uxx - Uyy)ux 



>'xy\, 



(5) 



where Tj 2 (^^^^2/ ~^ Vxxyx) 2^^^yyy ~^ ^yxyy) 

^ijlxxxy 'Hxxyx^ — 2^'^'^yyy "^yxyy)' -^^ &lS0 followS 



M 



that, up to boundary terms, (5) is equivalent to: 



5S, 



H — j d xdt^xxxyi^^xx ^yy^'^xy 

H" 'Hxxxyi^xx H~ '^yy)^xy\' (6) 



3 



To obtain this, we have used the faet that, up to bound- 
ary terms, rj^xyx does not contribute to the action. 



Phonon Hall viscosity as a response property of 
the electron system 



to Berry curvature^"^ of the ground state wave function. 
For a Hamiltonian ff[{Ai}] that depends on a set of pa- 
rameters {A,}, we have 



dH \ _ dE_ 



(11) 



In the adiabatic approximation, it is assumed that the 
motion of the lattice is infinitely slow compared with the 
motion of the electrons, so that at any moment, the elec- 
trons are in their ground state with respect to that par- 
ticular instantaneous configuration of the lattice. Within 
this approximation, the effect of lattice displacements 
is to alter the parameters in the effective tight-binding 
model for the electron system. The dependence of these 
parameters on the atom positions can be calculated us- 
ing standard ah initio methods. Thus, in the adiabatic 
approximation, the electrons will be described by an ef- 
fective tight-binding Hamiltonian, where the lattice dis- 
placements appear as external parameters: -fft.h. [{^i}]. 
Un is the displacement of the nth atom from its original 
location. We may also view Ht.b. as a function of the 
Fourier components Uq = -;^= u„e'^'". Then, the 
following linear response formula 
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die'' 



duo 



du- 



q,b 



(7) 



gives an additional term to the acoustic phonon effective 
action of the form 

^S=\j d^+^xd'^+^x'vabix - x')u^{x)ub{x'), (8) 

where x here is a (d + l)-component vector including 
space and time, r]ab{x) is the Fourier transform into real 
space-time of rjabiq, and we have taken the continuum 
limit to get the displacement field u{x). The leading 
order term that is independent of uniform displacements 
w ^ li -|- a is given by (2). Starting from (8), we find 



riijkl = \ lim lim 



d 



(9) 



For spatially homogeneous deformations, the distortion 
tensor Wij = diUj is a constant. To calculate riijki, it will 
be more convenient to take Wij to be a constant and to 
treat it as a parameter in Ht.b.- Then, 



Vijkl =1 lim 



cold J 



dte 



,iuit 



dH,.b.^^^d_H^^^^ 



dw 



'kl 



(10) 



For spatially inhomogenous deformations, we can con- 
tinue to use the DC response (10) instead of the exact 
AC response, as long as the acoustic phonon frequency 
is much less than the electronic energy gap. 

It is well-known that the adiabatic response of a Hamil- 
tonian to changes in some parameter is directly related 



where Ojj is the Berry curvature of the ground state wave 
function. Thus rjijki is given by the Berry curvature as- 
sociated with adiabatically varying the distortion tensor 
Wij = diUj as external parameters: 



1 



d 



^ l/\ ^ 



d 



dwki dwi 



■IV') + {i<^k) 
(12) 



where IV') is the ground state of the tight-binding Hamil- 
tonian Ht.b.- 



III. EXAMPLES 

In this section, we will study some examples of systems 
with a phonon Hall viscosity, including electrons hopping 
among the s-orbitals of a square lattice in a background 
magnetic field, a simple model for the quantum anoma- 
lous Hall state in HgMnTe quantum wells, and a simple 
mean-field model of a spinless px + ipy superconductor. 
In ccirtain limits we compare the phonon Hall viscosity 
of these systems with their conventional Hall viscosity 
studied in the literature. 



A. Hofstadter Model 

Consider a square lattice with nearest and next-nearest 
neighbor hopping: 



(13) 



Consider hopping among s-wave orbitals, in which case 
tij and iij depend only on the distance |rj— r^j be- 
tween atoms. To leading order in the crystal defor- 



mations, ti 



i-t-x 



t ~\~ t llxxj ^«,z-t-y — t ~\~ t Uy 



and 



UA+±±y = t + t'{j{uxx + Uyy) ± Uxy)- If t{r) is thc hop- 
ping matrix element between the s-wave orbitals that are 
a distance r apart, t = t{a), i = t{\/2a), t' = a^\a, 
i' = v^af^ly^a- In the absence of a backgromid elec- 
tromagnetic field and for constant lattice displacements, 
the Hamiltonian can be written in momentum space as 
H = X]k^kcj(.Ck. We note that a significant effect of thc 
strain fields on the energy of the system is to change the 
on-site energy of atomic orbitals. However, this contribu- 
tion does not affect the Hall viscosity below, so we ignore 
it. 
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In the continuum limit, the dispersion is, up to a con- 
stant, 

Cfe ^ -^hkjQij - {f + i'){Ua:x + Uyy), (14) 

where the effective mass is defined by = {t/2 + t), 
and Qij = 6ij + Sgij, 



6g = 2m 



I 2Uxv 



2ux 



(1 + J7)Uyy + Ux 



(15) 



In the presence of a gauge field, wc take k —iD 
—i{d — iA), so the effective theory becomes 



H = 



1 



2m 



-QijDiDj - {t' + i'){UxX + Uyy). (16) 



The phonon Hall viscosity is then related directly to the 
gravitational Hall viscosity of the electronic fluid^*. We 
find 



ss 

ixxxy 



it/2 + t)-^^-rjf^. (17) 



For Nl filled Landau levels, ry^ = Nitin/A, where n is 
the density of electrons, is the Hall viscosity of the elec- 
tron liquid when the crystal is ignored. ^"^'^^ The prefac- 
tor {t/2 + i)~'^i't' can exphcitly be verified to be of order 
one for typical s-wave orbitals and typical separations 
between atoms; it can also be calculated fairly precisely 
using ab initio methods. Note that since the Hamiltonian 
depends only on the strain field Uij , the phonon effective 
action also only depends on Uij and there is no depen- 
dence on the rotation tensor mij; i.e. rj^^ = r]^^ = 0. 

Therefore, we see that the result of Ref. 14 and 24 for 
integer quantum Hall (IQH) states has a direct effect in 
the phonon response, which, as will be discussed in Sec. 
IV, is a directly measurable physical quantity. For higher 
density systems, we cannot take the continuum limit; the 
phonon Hall viscosity is still a well-defined quantity that 
is calculable through linear response theory, but the pre- 
viously defined "gravitational Hall viscosity" of Ref. 14 
and 24 is not well-defined. 



B. Quantized anomalous Hall state and 
Hgi yMnj,Te quantum wells 

Here wc will calculate the phonon Hall viscosity for 
Hgi_j^Mnj^Te quantum wells, which exhibit a quantized 
anomalous Hall (QAH) state for certain thicknesses of 
the quantum well and spin polarization of the Mn ions^^. 
The QAH state is a band insulator that exhibits a quan- 
tized Hall conductance in the absence of a net magnetic 
field. The first lattice model for such a state was intro- 
duced by Haldane^^, and since then it has been proposed 
to be realized in Hgi_yMnj,Te quantum wells. ^^"^^ As the 
quantum well thickness and the magnetization of the Mn 



ions is tuned, the system can be tuned between differ- 
ent topological states: a quantum spin Hall state, QAH 
states, and the topologically trivial state. 

At the topological phase transitions, the phonon Hall 
viscosity exhibits non-analyticities that can be accounted 
for in the continuum Dirac approximation. In what fol- 
lows, we will calculate the phonon Hall viscosity for phys- 
ically realistic parameters, we will isolate the universal 
contributions that depend only on the low energy physics 
near the Dirac cones and we make contact with the cal- 
culations of Ref. 16 for the regularized gravitational Hall 
viscosity of the continuum Dirac model. 

The model for Hgi-yMn^Te quantum wells is given by 
a four-band Bloch Hamiltonian: 



H{k) 



h+{k) 
h-{k) 



(18) 



where the two-band Bloch Hamiltonians can be expanded 
in terms of Pauli matrices h±{k) = e±(fe)I -|- d±{k) ■ a, 
and h-{k) = h*^{—k). In the continuum limit and in the 
absence of lattice distortions, expanding near the F point 
k = (0, 0), we have 



yh 

2 , ;,2 



d±,x + id±^y = A{±kx + ik 

d±,.=M±-B{ki + kl„ 
e{k) = C±- D{kl + kl), 



(19) 



where M± = M ± 5M/2. The parameters A, B, C, D, 
M, and the lattice spacing a are given in Ref. 30 and 31 
for HgCdTc/HgTc quantum wells and the relevant ones 
are listed in Table I. 5M depends on Mn doping and spin 
polarization, as discussed in Ref. 25. 

The phonon Hall viscosity will be a sum of the contri- 
butions of each of the two blocks: 



Vijkl — Vtjkl + ''lijkl ' 



(20) 



where 



=2 8^^ I 



dd± 

d{diUj) d{dkUi) 



+ {i^k) 



(21) 



In order to calculate rjijki , we need to obtain the Hamil- 
tonian as a function of lattice distortions. To do this, 
observe that the blocks h± (fe) are composed of the spin- 
orbit coupled states |,s,±i) and Ip^jiipji; ±|)-^° Concen- 
trating on a single 2x2 block - for definiteness consider 
/i+ - the Hamiltonian is written as: 



(22) 



where i = x,y and n labels the sites of a two-dimensional 
square lattice. <r is the vector of Pauli matrices and the 
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hopping parameters are, to first order in lattice distor- two bands is: 
tions, 

ti = t + at'diUi, 
ti = t + at'diUi, 

= i{X + aX'dxU.^){x + d^Uyij) 
By = i{X + aX'dyUy){dyUxX + y) (23) 



The lattice parameters A, t, t, and m used above are 
related to the continuum parameters A, B, C , D, and M 
through: A = A/ a, M± ^ m± + 2t, B = C = 2i, 

and D = t/2. The hopping parameters are functions 
of the distance between neighboring atoms; the prime 
indicates a derivative with respect to this distance. The 
upper 2x2 block has topological phase transitions as 




2B 

otherwise 



(24) 



and similarly for the lower 2x2 block. 

From the Kubo formula (21), we see that the only non- 
zero terms are rj^^^y = —rj^yy^- The full effective action 
is given by (6); making the physically reasonable approx- 
imation aA' ~ A and at' « t, we find 



Vxxxy 



h r f aA a'^M+ 

~J ' - — • - — 



Sna^' \2B' 2B 



(25) 



2B 



is tuned. The Chern number Ci of the lower of the where / is a function of dimensionless parameters: 
I 



/(a,/3) 



sin^ kx:{P — 2 + cos ky)(fk 



[Q;2(sin^ + sin^ ky) + {fi - 2 + (cos k^ + cos fcy))2]3/2' 



(26) 



In Fig. 1, we plot the function -^/{a, /3) as a function 
of /3 for various choices of a. Note / has non-analyticities 
at the quantum phase transitions /? = 0, 2, and 4. 

The full Hall viscosity is given by the sum of contri- 
butions from the two blocks (see 20), which corresponds 
to taking the difference of the curve in Fig. 1 for two 
different values of m: 



(27) 

as a 
m+ — 

m_. In Figs. 2 - 3, we plot rj^xxi, as a function of either 



where f± 



f ( Aa a^M± 
■I \2B' 



2^ ,. We may view rj^^^y 
function of M = {m++m-)/2- 4B/a'^ and SM 

']xxxy 




FIG. 1. The function -^^ 
for various choices of a. 



rf{a,l3), plotted as a function of , 



A (eV ■ A) B (eV ■ A^) a (A) M (eV) 



3.645 



-68.6 



6.46 -0.01 



TABLE I. Realistic parameters for HgTe quantum wells, 
taken from Ref. 31 



M or SM, using physically realistic parameters (Table 
I) and focusing on the region near one of the transitions 
into the QAH state. 

In Fig. 2, we fix M and plot rjxxxyiSM); experimen- 
tally this can be done by tuning an external magnetic 
field. We observe a discontinuity in the slope of rj^xxy as 
a function of 5M at the transition into the QAH state. 
We denote this discontinuity A ^g^^f " • ^^S- ^' 
a?5M/2B <C 1 and plot rj^xxyiM). This shows a disconti- 
nuity Ari{M). While a full lattice calculation is required 
to compute the Hall viscosity, the properties of these non- 
analyticities can be accounted for in the continuum Dirac 
approximation to the above lattice model. In this approx- 
imation, we take sinfc ^ k and cosfc « 1 — k'^/2 in (26); 
near the topological phase transition at /3 = 0, we have; 



27ra^ 



(e-1) 



(a2 - e)2 



0(A), 



(28) 



where A is a high-energy cutoff. While the Hall viscosity 
will in general depend on A, the first term above is re- 
sponsible for the non-analyticities of / and consequently 
of Vxxxy Using (28), we can estimate the discontinuity 
in drjxxxy/9SAI at the transition, as shown in Fig. 2: 



(^Vxxxy 



dSM 



h \B\ 



TT a'^A'^ ' 



(29) 
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Similarly, the discontinuity in r]xxxy{M) for fixed 
a?5M/2B <C 1 at the transition is found to be 



h B 



6M 



(30) 



Finally, we note that the discontinuity in the derivative 



of [M) can also be computed using (28): A 



dM 



OC 



[q^2 |e=0+ 



A 



xxxy 



dM 



We find: 

ft . 
i -5M 



2 

l2 



(31) 



The first term above depends only on parameters of the 
low-energy theory of the state and is independent of the 
lattice spacing and other high-energy details. This is es- 
sentially the contribution that was found for the gravita- 
tional Hall viscosity in the regularized Dirac model stud- 
ied in Ref. 16. We see that the other contributions to the 
Hall viscosity that we find are dependent on the lattice 
spacing and other high energy details, which is why they 
are missed in the regularization of the continuum Dirac 
model of Ref. 16. 



C. Interacting states and + ipy superconductors 

For interacting systems, while there is no universal re- 
lationship between phonon Hall viscosity and the gravita- 
tional Hall viscosity, a major exception occurs in systems 



with only on-site interactions: H ~ J^ai'^' 



'c]ci 



h.c.) + U^^rii^ni]^. If the hopping involves s-wave or- 
bitals, then in the dilute limit, where the system can be 
described by a continuum interacting theory, the effect 
of a strain in the lattice is equivalent to a deformation of 
the gravitational metric. For such systems, the phonon 
Hall viscosity is then directly related to the gravitational 



' (A 




FIG. 3. Plot of ■na:xxy/hSM as a function of M for fixed SM = 
W~^eV and for realistic parameters (Table I). There is a 
discontinuity at the transition to the quantum anomalous Hall 
state, which occurs at SM = ±M. 



Hall viscosity, through a proportionality factor of order 
unity, as in (17). 

As an example of an interacting state with a phonon 
Hall viscosity, we consider a BCS mean-field description 
for a px +ipy superconductor, since such a model may be 
relevant for the chiral superconductor Sr2Ru04^2 First 
consider the interaction between the nearest and the next 
nearest neighbors, 

U^~VJ2 4c.4c, - F ^ cUc]c,. (32) 



Taking only the Cooper channel of the interaction yields: 
= - Vkk'4/clk'C-kCk, (33) 



k,k' 



where 



kk' 



N 



[2{v + v) + (r + r){uxx + uyy) 



{V/2 + V)(h-k[){k,-k'^yg,,, 



(34) 



200 300 

6M (meV) 



FIG. 2. Plot of rjxxxy/h as a function of 5M for fixed M and 
for realistic parameters (Table I). There is a discontinuity in 
the slope at the transition to the quantum anomalous Hall 
state, which occurs at 5M = ±2M. 



gij = + 2{v/2+v) ^^ij has the same form as g^j in eq. 
(15), with t and t replaced by V and V: 



\l + ^)Ux 

2u 



2Ur 



xy 



J.yy ^U,xy 

v 



(1 + T77)"to 



(35) 



V = a^\,=a and V' = V2a^|^^^^, where V{r) and 
V{r) are the nearest and next-nearest neighbor interac- 
tions, which only depend on the distance r between the 
nearest or next-nearest neighbor sites. As will be dis- 
cussed below, the phonon Hall viscosity is proportional 
to the gravitational Hall viscosity only when the two met- 
rics, gij and gij are the same. 
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To study the simplest possible scenario, we assume the 
electrons hop among a single s-wave orbital of the atoms. 
The BCS mean-field Hamiltonian is then 

iJscs = ^(ek-/x)cj,Ck + i^[Ak4clk + ^f.c.], (36) 

k k 

where €k is given by Eq.(14). For + ipy pairing, we 
take the order parameter to be, 

Ak =A(sin kx + i sin ky) 

+ Aisin kx cos ky + i sin ky cos kx) ■ (37) 

The order parameter must satisfy a self-consistency equa- 
tion: 



Al 



Ak' 



r 

k 



(38) 



where Vkk' is the Cooper channel of the interaction and 
Ek = V(ek - m)2 + |AkP. 

To calculate the phonon Hall viscosity for this system, 
we need to obtain the effect of the lattice deformation on 
the order parameter A/j. For simplicity, we consider the 
long- wavelength continuum limit, where 



(39) 



where = (4,c_fe), and 



BdG 



A*{\k\){kiex^ - 



ikjCyj) 



A{\k\){kiexi + ikjCyj) 
-ek + iJ- 

(40) 

ki is a unit vector, and A(|fc|) is chosen to fall to zero 
far away from the Fermi surface (the cutoff scheme is 
explain in Fig. 4), /i is their chemical potential, and e^j = 
5ij + Scij where Scij is a linear combination of the lattice 
distortions wu = dkUi. It is convenient to define an 
"order parameter metric" = CiaCja, which we can 
fix in terms of g and g using the BCS self-consistency 
equation. In the continuum limit, where the system is 
rotationally invariant, it is simple to show (see Appendix 
A) that 



19ij + (1 - l)9ij, 



(41) 



where 7 is a constant that can be determined from the 
self-consistency equations (see Appendix A). To first or- 
der in Wij, we have gf^ = Sij + {Scij + Seji), so the above 



Scji 



equation does not fix Scij — Seji . We may fix Se^ 
by observing that the only affect of a rigid rotation of 
the crystal should be to rotate k^ and ky into each other. 
Thus: Cij — eji = 2my. These considerations fix the 
dependence of the order parameter on the lattice defor- 
mations. Thus we can now use the Kubo formula and 
explicitly obtain the Hall viscosity: 



Vijkl 



1 h f 

28^7 



(fkd ■ ( 



dd 



dd 



dwij dwki 



) + {i^k) (42) 



In Appendix B, we present some details of the calcula- 
tion. In the case where gij = gij, the calculation simpli- 
fies considerably and we find a simple result 



i'l' 



1 



»7 



-hn 



ft' 



r\9 'igr'> 



(43) 



where the gravitational Hall viscosity is ry^ = jhn. The 
constant of proportionality rj^ /rj^r typically of order 
one. In Fig. 4, we show the Hall viscosity calculated 
from the results presented in Appendix B for gij 7^ gij. 
While this does not modify the Hall viscosity in the weak 
pairing limit much, the behavior close to the weak to 
strong pairing transition is dependent on completely non- 
universal features, such as the frequency dependence of 
the pairing gap. 

For a system with a circular Fermi surface, the Px + 
ipy state has another interesting feature: the U{1)^ x 
U{1)l^ symmetry associated with particle number and 
angular momentum conservation is spontaneously broken 
to a diagonal subgroup, ?7(1)l^_0. This implies that in 
the effective action, nixy and 0, the angle through which 
the crystal is rotated about the z direction and the overall 
phase of the order parameter, respectively, should appear 



together as 



-(p. This means that the effective action 



of the crystal involves the phase of the order parameter 
as well: 

SS,=2jd^.dt[,-iUxx-Uyy)Uxy 

+ V {"^xx + ^yy){^xy 

Physically the U{1)l^-4, symmetry requires to in- 
crease by h when two electrons are adiabatically added. 
Since m^y and (f) are conjugate variables to and the 
number of Cooper pairs, the above effective action con- 
tains the Berry phase term for this adiabatic process: 



dL 



H 



drrix 



= SL,. 



(45) 



The physical consequences of the coupling between the 



uniform compression, 



and d) means that the 



change in particle density is given by 



Sn = 4?7 (uxx + Uyy) = 4r? 



M 



SA 



(46) 



where A is the area of the 2D system and n is the particle 
density, when n is held constant^^. A generic supercon- 
ductor of course has 



Sn/n = a{fi)6A/A, 



(47) 



where a{ij,) is a constant depending on the chemical po- 
tential; note that we would have a = — 1 if the total par- 
ticle number had been fixed. The additional symmetry in 
this problem, which relates phase rotations to rotations 
of the crystal, sets a(/x) = —Aij'^/n. This suggests a 
way to experimentally measure r]'^ in systems that have 
an additional symmetry involving spatial rotations of the 
crystal. 
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FIG. 4. The px + ipy superconductor phonon Hall viscosity. 
The filled curves are for the cases with gtj ^ gtj, while the dot- 
ted curves are for gij = gij. We have set V' /V = 2t' /t, t/t — 
-t'/t = 0.75, and V/V = ~V' /V = 0.40. The pairing gap 
was set to have a Gaussian cutoff A = Afcexp[— fc'^/2(5fc)^], 
with Sk = 1/a and A(2me//a/fi,^) = 0.3,0.4,0.6 for the 
blue, purple and red curves respectively. Note that while 
rj^ converges to the same asymptotic values in the weak- 
pairing limit, in the strong-pairing phase and near the quan- 
tum phase transition, it is not proportional to the density n 
when gij ^ gij. 



IV. PHYSICAL CONSEQEUNCES OF PHONON 
HALL VISCOSITY 



A. Acoustic phonon dynamics 

Consider the effective long-wavelength elasticity the- 
ory of a crystal, given by (1) and (2). Note that this 
is an expansion in the displacement fields and its gra- 
dients. Since for the sound waves uj oc + • • • , the 
Hall viscosity terms are actually of order k^, so for 
consistency one must also include a term of the form 
6S3 = J d'^xdtXijkimdmdiUjdkUi, but such a term van- 
ishes in the presence of inversion symmetry. As we noted 
in Section II A, anharmonic effects are 0{k'^), so Hall 
viscosity may be distinctly measurable because its effects 
appear at lower order in k. We briefly mention the effects 
of impurities later. The physical consequences of Hall 
viscosity terms can be analyzed most simply by consid- 
ering 2D systems whose long- wavelength elastic theory is 
isotropic. This would be directly physically relevant for 
2D systems with square lattice symmetry; the considera- 
tions directly apply for layered 3D crystals as well, where 
the 2D layers have a square lattice symmetry and where 
we consider phonons with wave-vector oriented parallel 
to the 2D layers. For such systems, the elastic theory 
simplifies and one obtains for the equation of motion: 

ii, = c^V^u, + {cf ~ c^)d^V ■ u -I- rf^'^e'^iij/p, (48) 

where the indices i, j run over the 2D spatial coordinates, 
Ct and ci are the transverse and longitudinal sound veloci- 



ties, respectively, and 77 = rj^xxy is the Hall viscosity. It is 
simple to show that such a wave equation does not admit 
purely transverse or purely longitudinal solutions. Let us 
denote e± as the eigenmodes of the system, and let the 

basis ^1 0^ and ^0 1^ correspond to the longitudinal 

and transverse acoustic phonon modes, respectively. In 
this basis, the eigenmodes of the system in the presence 
of the Hall viscosity are 



e+ cx 



where x 



1 




+ (49) 



This defines the characteristic frequency 



is a dimensionless parameter. 

n 

where B is the bulk modulus of the crystal. (Note 
that for a crystal, q > acn for some constant a of order 
unity, which is why p{c^ — cf) ^ pcf ~ B). Observe that 
to hnear order, there is a ^sgn{r]) phase shift between 
the longitudinal and transverse modes. The dispersion 
relation is 

2 ! 2 , 2 , 2 ; 2 / 2 I 

oj =y[q +Ct +f] k /p ± 

yjc^ + (c2 + k^ri-^/p-^Y + 2c2(fc2r;2/p2 _ ^2)]. (59) 

The shift in frequency for a given acoustic phonon mode 
for finite 77 is Aui/iu{r] = 0) ~ x{uj{r] = O))^. In principle 
then the shift in frequency can determine r/. However, 
since this is not sensitive to the sign of the Hall viscosity, 
it may not be a useful method in practice for determining 
the Hall viscosity. 

An analysis of surface (Rayleigh) waves of a 3D 
medium with non-zero rjxxxy = Vyyyx displays similar be- 
havior. For a medium with surface at z = 0, a surface 
wave travelling in the x-direction must have Uy — in 
the absence of Hall viscosity, due to stress-free boundary 
conditions. In the presence of a Hall viscosity, the surface 
wave acquires a Uy component, which to linear order in 
rjuj/ pc^ differs by a phase shift of 7r/2 and has a relative 
amplitude of rjio/ pc^. 



B. Numerical estimates and discussion of possible 
experimental detection 

As explained above, the physical consequences of a Hall 
viscosity in the phonon effective action is the mixing of 
longitudinal and transverse sound modes, which at a fre- 
quency ijj is determined by ijj/ujy, with LOy = pc? jrp a 
characteristic frequency scale associated with the Hall 
viscosity. To lowest order in u /lo^^ there is a phase shift 
of 7r/2. For the bulk modes, a more precise value of uj^ 
is LOy — p{cf — cf)/r]^ , and the amplitude ratio between 
the two modes in the elliptical polarization is uj/tUy. It 
is not clear what the minimal experimentally measurable 
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value of the Hall viscosity is. From elementary consid- 
erations, we can put a rough bound on what may real- 
istically be measured. First, for the sound waves to not 
destroy the crystal, we expect that the strain is small: 
\du\ <C 1, which implies that \ku\ <C 1. Since we have 
roughly a; ^ cfc, this implies wu ^ c. The amount of mix- 
ing is determined by ui/cOy, so the amount of amplitude 
from the other mode that is mixed in is {uj/uJv)u. For 
this to be realistically measurable, this should be much 
larger than the size of the quantum fluctuations of the 
wave function of an atom, which is on the order of 0.1 
A. That is, \lliu\ ^ ujy x O.lA. Thus we have the condi- 
tions uiy X 0.1 A <^ u <^ c. For a typical sound velocity of 
5 X 10^ cm/s, this implies 

a;„ < 5 X 10"s-^ (51) 

For smaller sound velocities, this bound will be smaller. 
This is not a fundamental bound, but a practical one. 
This is because in principle it is possible to measure os- 
cillations of the center of the wave function of an atom 
at a resolution that is smaller than the characteristic size 
of its wave function. 

The two-dimensional mass density of the crystal is p a 
Amp/a^, where A is the atomic number of atoms of the 
crystal. Typically, A ~ 10, and o ~ 4 x 10~^cm, so p ~ 
10~^g/cm'^. Furthermore, cf ~ ^ {a^ — 1)10^"cto^/s^, 
where a = ci/ct and typically a ~ 2. r;^ ~ rj^^ ~ hUg, 
where rig is the electron density. Thus, for a typical 
2D electron density of lO^^cm"^, with a ~ 2, we see 
iJv ^ lO^^s"^. For a 1 GHz measurement, oj/cOv ^ 10~^; 
at 100 GHz, io/uiv ~ 10~^ and, depending on material 
parameters, could be closer to 10^^. Note that the fre- 
quencies must be much less than the energy gap of the 
electronic state, which for a 10 K gap translates to ap- 
proximately 2 THz, and also much less than the phonon 
Debye frequency, which is close to 10 THz. 

For quantum Hall states induced by an external mag- 
netic field, the electron densities are usually low, rig ~ 
lO^^cm"^, yielding immeasurably small values for co/u)y. 
An exception may be graphcne, where recent advance- 
ments in applying extremely large gate voltages may 
allow for much larger densities of electrons participat- 
ing in quantum Hall states.'^* The necessary values of 
Ue ~ lO^^cm"^ usually appear in states that sponta- 
neously break time-reversal symmetry, such as quantum 
anomalous Hall states, ferromagnetic insulators, or chiral 
superconductors, where the effective magnetic moment 
per lattice site is much larger than could be produced by 
an external magnetic field. 

The effects discussed here would most easily be mea- 
surable in bulk, layered 3D crystals, for phonons propa- 
gating along an in-plane high-symmetry direction. While 
there are a number of examples of 3D IQH states, the 
value of the Hall viscosity is probably too small to be 
measured, since the particle density is too small, though 
not typically as low as in 2D quantum wells. More 
promising systems are those that spontaneously break 
time-reversal symmetry, because those typically will have 



much higher angular momentum densities. One promis- 
ing candidate may be the chiral superconductor Stron- 
tium Ruthenate, which exists as a 3D crystal and may 
have a large enough Hall viscosity because it sponta- 
neously breaks time-reversal. Another promising set of 
materials to measure a phonon Hall viscosity are 3D 
ferromagnetic insulators, for example those discussed in 
Ref. 35-37 . Note that in cases where the spin gap is 
small, an external magnetic field can be used to ensure 
the electronic state is fully gapped. 

In principle, one way to measure such an effect would 
be through pulsed echo ultrasound mcaurements, which 
have been successful in detecting circular polarization 
between transverse sound waves. However, while bulk 
pTilsed ultrasound seems to be limited to frequencies on 
the order of 1 GHz, it is not clear what the ultimate 
bounds are on an experimentally accessible amplitude ra- 
tio between transverse and longitudinal waves. A more 
promising experimental technique appears to be time- 
dependent x-ray diffraction'^* . Such techniques have been 
developed only recently over the last decade and have 
been used to directly image acoustic phonon modes^^'^°. 

One complication of measuring the phonon Hall vis- 
cosity is related to the effects of crystal disorder, which 
can also mix transverse and longitudinal waves. How- 
ever, the effects of disorder are not sensitive to the sign 
of the time-reversal symmetry breaking of the electronic 
state; this dependence on the sign of the time-reversal 
symmetry breaking is unique to the Hall viscosity, and 
can be used to extract the phonon Hall viscosity even for 
imperfect crystals. 

We would like to point out that related phenomena 
occur in various other time-reversal breaking systems. 
The phonons in a ferromagnet, for example, can ex- 
hibit acoustic Faraday rotation, where the two trans- 
verse modes acquire a circular or elliptic polarization^^'^^. 
However such systems cannot be described by a simple 
local effective action in terms of the strain fields because 
they are coupled to magnons, which are gapless; inte- 
grating out the magnons will result in non-local terms in 
the crystal effective action. The physical manifestations 
of such phenomena are also quite different; they occur 
as resonances when the frequency and wavelength of the 
phonons and magnons are matched. More directly re- 
lated phenomena have been considered in the case of ionic 
crystals in an external magnetic field^^, and in Tkachenko 
modes of vortex lattices in rotating superfluids*^"^®. In 
these situations, one obtains a related equation of motion 
as in (48). 



V. CONCLUSION 

We have proposed the acoustic phonon Hall viscosity 
as a novel probe into the adiabatic Berry curvature of 
the many-body electron wave function for gapped, time- 
reversal symmetry breaking electronic states, and we 
have computed it and studied its behavior in a number of 
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theoretical models. In some simple cases, we have found 
that the phonon Hall viscosity is proportional, with a nu- 
merical factor of order 1, to the gravitational Hall viscos- 
ity of the continuum electronic theory. Our numerical es- 
timates indicate that this is a measurable effect and there 
may be a number of materials, particularly the ferromag- 
netic insulators, which might be suitable candidates for 
experimentally detecting the phonon Hall viscosity by 
measuring time-reversal symmetry breaking corrections 
to acoustic phonon dynamics. Since the effects are ex- 
pected to be small and their measurement would require 
high spatial resolution, it appears that time-dependent 
x-ray diffraction may be currently the most promising 
probe. As phonon Hall viscosity is developed into a more 
mature experimental probe, we hope that it can eventu- 
ally be useful as a novel lens into the possible topological 
behavior of electron systems. 
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Appendix A: p + ip BCS self-consistency 

The BCS mean-field Hamiltonian is: 

HBdG = "^{^k - M)CkCk - ^ AkC_kCk, 
k k 

where 

Ak =A(sin kx + i sin ky) 

+ A (sin kx cos ky + i sin ky cos k^) 



(A2) 



and the kinetic energy is that of Eq.(14). Ak satisfies a 
self-consistency equation: 



Ak = - ^ li.k' 



Al 



(A3) 



where Vkk' is the Cooper channel of the interaction and 
^k = V(ek-/0' + |AkP. 

When the crystal is strained, the order parameter will 
take the following form in the continuum limit: 



A/j — ^{kiCxi 4~ ikiCyi)^ 



(A4) 



where eablidiUj}] are functions of the distortion tensor 
Wij = diUj. In order to obtain the phonon response, we 
need to obtain this function to linear order in diUj. First, 
observe that due to the U(1)l^ x U(1)0 — >■ U(1)l^_0 sym- 
metry breaking, (j) spatial rotation is equivalent to the 
gauge transformation A Aexp{i^). Therefore 



Next, observe that e^x, ^yyj and e 



^yx' 



(A5) 

xy -\- Byx are symmetric 
in X and y, so they can only depend on the strain tensor 
Uij. The self-consistency equation can be thought of as 
a constraint on the "order parameter metric" g^: 



Considering the variations of this: 



(A6) 



df 



5g 



dj_ 
dgij 



Of 



dgi. 



(A7) 



oc 5ij, and 



For a rotationally invariant system, df /dgij 
similarly for df /dgij. This implies that 5g^ = "fSga + 
^Stjii, where 7 and 7 are constants. Now observe that 
when gij = gij, we are merely implementing a coordi- 
nate transformation, so we should have gfj = gij, which 
implies 7 + 7 = I. Furthermore, for a rotationally in- 
(Al) variant system, the deformations dgxx = —Sgyy = e and 
Sgxx = —^9yy = e are equivalent to the deformations 
Sgxy = e and 6gxy = e, because for a rotationally invari- 
ant system, the two types of deformations simply differ 
by a rotation. Thus, we conclude: 



9^ = 19i3 + (1 - l)9ij- 



(A8) 



The constant 7 can be found from the self-consistency 
equation. 

To actually calculate 7 from the self-consistency equa- 
tion, we note that the assumption we made above tells us 
that the effect of change in the kinetic metric gij should 
be proportional to the effect of coordinate transforma- 
tion 5ij — >■ gij. This means that, if we consider the case 
^gxx = —59yy = ei and 5g = 0, the change in the order 
parameter should come out as ^Ak = 7eiA(/Ca; — iky)/2. 
Thus, to the self-consistency condition 



gAk 

de\ 



I 5Ak' Ak' / Ck' — M de 



+ 



d\A^ 



2Ek' dei 2El, \ E^' dei E^' dei 



(A9) 
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we can insert 



dei 



(AlO) 



to obtain 



Ak'|A'k|5JAk^ 
9ei 



(All) 



and solve for 7. For this step, it is convenient to eliminate d\Aw\/dei through a coordinate transformation argument, 
which makes use of the fact that in the integral over k', we can do a coordinate transformation {k'^, k'y) — >■ {k'^,ky) + 
jei{k'^, —k'y)/2 without changing the value of the integral because the Jacobian of this transformation is one, to linear 
order in ei. This gives us 



AknA'i,|9|Ak 
2El, dei 



1 dAu> Ak' ek' - n dei, 



2Ei,, dei 2El, E^, dei 



avkk' Ak> 



k, -^k' k' ^^^"^ """k' k' 

We can obtain the derivatives on the right-hand side by noting that this coordinate transformation gives us 



(A12) 



Ak' ^ Ak' + -7A{k'^ - ik'y), ek' + jeu{kl - k^), Fkk' ^ Vik' - lV{k^k'^ - kyk'y) 



(A13) 



where we used Vkk' = — 2F(k • k'). Since the change in |Ak'| is same, we can insert this result from the coordinate 
transformation into Eq.(A9). 

I 



We find at the weak coupling limit 
1 - ln(2ee/|Ao| 



l + ln(2ec/|Ao| 



(A14) 



where Cc is the cutoff energy of the pairing and |Ao| is 
the value of A at the Fermi surface. The negative sign 
of 7 is because our kinetic metric deformation increases 
the density of state around (0, ±fc/) compared to (±fc/, 0) 
and therefore it is energetically advantageous to have the 
Py pairing to be stronger than the Px pairing. 



Appendix B: p+ip BCS Hall viscosity calculation 

Now we would like to calculate the Hall viscosity for 
the Px + ipy BdG Hamiltonian. We want: 



167r2 



d'^kd- 



dd dd 

X 



dWr 



dUxy 



From HsdC = d - we have 



2m 



(B2) 



note that in a flat metric, d\\ = yjdx + dy and dz would 

depend only on k. In this section, since d|| has explicit 
dependence only on Cij while dz has explicit dependence 
only on gij , we will treat eij and gij as independent vari- 
ables. 

For our calculation, we want to write the derivatives 
in terms of g^j and Cjj, which gives us 



9d|| _1 / ad|| ^ ~9d|| 
d^xx 2 \ Ogxx dcyi 



— — (akx-^TT^ 
2\ -"dhx 



bk 



^ dky 



ddz 

dWxx 
^^xy 

ddz 

dllxv 



where 



ddz , ddz 
=a- h 0- 



dgx 
~ ( dd 



yy 



=26 



\ dCxy 

ddz 



dg 

-- ddn 
b- 



1 f , ddz 

2 r^^a/t. 



bk, 



dcy 



b h 



ddu 



ddz 
'dky 
dd. 



dgx 



^ dk " 



'dkx 



(B3) 



a ■■ 
b- 



rrieff {t' + i'), 
meffi', 



a = 7K//(t' + i')] + (l-7) 



V' + V 
V/2 + V' 



b = 7(me//F) + (1 - 7) 



V' 
V/2 + V' 



(B4) 
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and we used 



and 



94 
den 



94 



1 /. dd, 

ad 



-hi 



dk-i 



G%4 — 5% . 



To see this, first observe: 



dwi- 

dd, 
dwi. 



kl 



E 

kl 



dwij dcki ' 

dgu ddz 
dwij dgu ' 



(B5) 



(B6) 



Assuming Wxy-Wyx = Cxy-eyx, all dgki/dwij,deki/dwij 
vanish except for 



dgx 



dg. 



vv 



dwx 
dgx 



dWy 

dg. 



vv 



dW; 



yy 



dwx 



dgx y dg: 

dW; 



xy 



xy 



dw. 



(B7) 



yx 



a =2 
6=2 



dcx 



dWxx 

de 

^^xx 



vv 



dw, 

ydCxy 

dw. 



de 

^IJCjjy^ 

dWyy' 

2 '^^vv 



xy 

=2^^ + l 



dw. 



dWxx 

= 2^ + 1 

OWyx 

2 ^^yx _ ^ 



xy 

Now we need to calculate 
h f ,0, I dd 



dw. 



(B8) 



yx 



167r2 

h 



fed- 



dwx 



dd 



1 



-32.^-"^/'^v^• 



1 



ad" / dd' 

dkx \ dky 



327r2 

n 



d^k—d- 
d^ 



dd' 

"dkx 

dd 



k — X (k — k — 
^ dk,, \ dk,, " dk. 



32^2 fc.d k^d l^g^^x 



+ 



327r2 



kyk^d- 



dd' ^ dd' 



(B9) 



where d' = (6d||/6, rf^) and d" = (ady/a, rfx). Then, 
using 



d 
dkx 

d 
dk,, 



--COS( 



:sin( 



d sin (p d 



dk k dcf)^ 
d cos(j) d 



dk 



k d(t>' 



(BIO) 



we obtain: 



=Cr J k^dk^ + J k^dk 



i »x 1 dd ddz 



Ca / k^dk 



Idd 



dd^. 



=87rCin- -C2 



kdk{l - d^) + C2 / k'dkd,{l - df) 



1 dd 
'ddk 



C-i / k^dk 



d dk dk 



where 



Cl/h: 



1287r 



{ab + 3ab- 4bb), 
(a6 — 56). 



dd^ 
^ dk 

(dl 
\ dk 



^ddk 



1 

Udk^ 
(Bll)' 



(B12) 



Now note that when the kinetic and interaction metrics are equal, we have a = a, b — b, which means Ci/h = 
b{a — 6)/327r and C2 = C3 = 0. Inserting these into Eq.(Bll) gives us 



r]" = SnCm = 6(a - b)hn/4 = {t/2 + i) 



ft' hn 

XT' 



(B13) 
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In the weak pairing limit, even when the kinetic and interaction metrics are different, we have rj^ oc n under some 
reasonable assumptions. First, note that 

y kdk{l - rf^) = y kdk{l - 4) = Ann, 

in this limit. Meanwhile, if we assume |Ak| oc k near the Fermi surface, we have 

At the quantum critical point, the ratio of the phonon Hall viscosity to the gravitational Hall viscosity becomes 
different from what we have for the weak-pairing limit. We can see this from 



y d^ dk dk J d3 \ J Jo £+2rn;A2y 2 ^ ^ 



where dii = Afc. We also have 



47rn < y fcd/;;(l - d^J = j kdk{l - d^){l + 4 + ^1) < 127rn. (B17) 

In the strong-pairing limit, we do not find asymptotic limit to the ratio of integrals to n that is independent of cutoff. 
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